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^^ . We address a class of models in which neutrinos, having a small mass, origi- 

rvg . nate the highest energy cosmic rays interacting with the relic cosmic neutrino 

o: 

0^ ■ background. Assuming lepton number symmetry and an enhanced neutrino 
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,^ . density in arbitrary size clusters (halos), we make an analytical calculation 
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Q . of the required neutrino fluxes. We show that the parameter space for these 

c/3 . models is heavily constrained by horizontal air shower searches. Marginal 

^ . room is left for models with exceptionally flat neutrino spectral indices, neu- 

'^j . trino masses in the 0.1 eV range and supercluster scale halos of order 50 Mpc 
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size. Our constraints do not apply to models with lepton number asymmetry. 
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The detection of ultra high energy cosmic rays (UHECR) above the Greisen-Zatsepin- 
Kuzmin (GZK) cut off |[T|] has stirred the research activity in cosmic acceleration mech- 
anisms. Above the GZK cutoff protons rapidly loose energy through photoproduction in 
the cosmic microwave background and therefore sources must be relatively nearby 0. The 
known objects in our "extragalactic neighborhood" , within few tens of Megaparsecs, have 
difficulties to accommodate stochastic acceleration mechanisms (most commonly invoked) 
due to dimensional arguments. Moreover such UHECR deviate little in the magnetic fields 
encountered over this length scale and no obvious astrophysical candidates are seen in the 
arrival direction of the few detected events. Several possibilities have been considered to 
explain these events [0,H and in particular that cosmic ray production arises through ultra 
high energy (UHE) neutrino interactions with the cosmic neutrino background 0,^. 

UHE neutrinos could come from cosmic distances and interact with the relic neutrinos in 
our halo. The final stable products of these interactions would be gamma rays and protons 
(besides secondary neutrinos) which would constitute the high energy end (above ~ 10^^ eV) 
of the cosmic ray spectrum. The idea is attractive because it avoids the constraint that source 
candidates must be at distances below ~ 50 Mpc, but it requires large fluxes of very high 
energy neutrinos (above ~ 10^^ eV) ||^,|^ without getting into the details concerning the 
UHE neutrino production mechanism. Models involving annihilation of topological defects 
IPJ^ and heavy relic decays [0 could for instance produce these neutrinos rather naturally. 
Bounds for some of these models have already been discussed in the literature based on 
neutrino and photon flux measurements ||,|ll|-|T3[. 

Only the resonance peak in the Z^ production interactions with the cosmic neutrino 
background can provide any significant secondary particle flux. In the massless neutrino case 
an energy El^'^ ~ 10^^ GeV is required to produce Z°'s at resonance since relic neutrinos have 
energies ~ 2 K ^ f .7 10"^ eV. This possibility would imply either neutrino fluxes exceeding 
current limits from Horizontal Air Showers, as will be shown below, or unnaturally fine tuned 
neutrino energy spectra. If neutrinos are massive, a possibility that is becoming increasingly 
more realistic in the light of recent results by Superkamiokande [|l4l, the necessary neutrino 



beam energy to produce Z^'s at resonance in the interactions with the rehc neutrinos becomes 
inversely proportional to the neutrino mass. Moreover background neutrinos would tend to 
accumulate in an extended halo as pointed out in Ref. P,|6|, p!5| , p!6| increasing their local 
density with respect to the cosmological value and the probability of nearby interactions. 

The idea has already been discussed by several authors. Using an incoming neutrino 
flux of spectral index 2, Waxman has discussed the models from general energy density 
arguments, using limits on local neutrino density because of Pauli exclusion, to conclude 
that a new class of models would have to be invoked to accelerate the neutrinos themselves 
and that the energy required is comparable to the total photon luminosity of the Universe [§ . 
The calculation depends on the spectral index and the energy cutoffs of the assumed neutrino 
spectrum. Yoshida et ai. have recently computed the particle spectra for several case studies 
after detailed propagation of all secondary products in the extragalactic magnetic fields and 
through the cosmic microwave background, and assuming clustering in supergalactic scales. 
These cases support the possibility that the produced UHECR, neutrinos and gamma rays 
are compatible with neutrino observations and bounds ||^. 

In this paper we further discuss this idea analyzing the model dependence on the assumed 
neutrino spectral index. We establish the neutrino fluxes firstly using energetic considera- 
tions similar to those of Ref. |^ and then analytically calculating the proton and photon 
secondary spectra. If large fluxes of high energy neutrinos exist, they should have been de- 
tected. Indeed by imposing that the neutrinos produce the observed UHECR one expects a 
much larger flux of neutrinos as pointed out in Ref. , since the probability for interacting 
in the neutrino halo is small. In a phenomenological approach we leave the neutrino spectral 
index and the local neutrino density enhancement in the halo as free parameters. We will 
show that for a large region in the two dimensional parameter space, the required neutrino 
flux is heavily constrained by existing data on horizontal showers. 

Assuming neutrinos are massive (of order 1 eV) we consider a neutrino of energy E^y. 
This neutrino could interact with an antineutrino from the cosmic background, with a center 
of mass squared energy s = 2m,yE^. The cross section for this process is maximal near the 



Z^ resonance, of width Y z-, which occurs for neutrino energies of £"1, ~ 4 10^^(1 eY /mi) 
GeV. In hadronic decays the Z^ produces high energy particles, mostly pions, which further 
decay so that only high energy photons and protons (and neutrinos) would eventually reach 
the Earth. The final particle spectra are given by a convolution of the quark fragmentation 
functions and the in flight decays of all the intermediate particles. These spectra have to 
be propagated in the photon cosmic background, IR background, galactic fields, etc. which 
would alter the arrival fluxes of high energy photons and protons. As shown by Yoshida et 
al. 0, the final particle spectra agree well with observations and can explain the UHECR 
spectrum for a wide range of spectral indices, 7 ~ 2, assumed in the original neutrino flux. 
The survival probability of a UHE neutrino in the relic neutrino background is in general 
given by Vs{Ey) = e""^", with r^ being the opacity. Considering only the resonant Z° 



production cross section a^p [17| in a matter dominated Universe the opacity can be well 



approximated by the redshift integral: 

c /•^'"- dz{l + zY a^p[2mpEMil + z)] 

Tu (E^o) ~ -- ripo / , (1) 
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Here H is the Hubble constant, E^, is the interacting neutrino energy and the subscript is 
used to indicate the present value of a redshift varying quantity. Qm and ^a are respectively 
the matter density and the cosmological constant terms in the Friedmann equation expressed 
in dimensionless units. The energy integral over the relic neutrino spectral density has 
been eliminated in the assumption the neutrinos are non relativistic so that rip is the relic 
antineutrino number density and the argument of the interaction cross section is the redshift 
varying center of mass energy of the collision. 

If one integrates this expression to galaxy formation era Zmax ~ 5 assuming no cluster- 
ing of the relic neutrinos, the uncertainty in the numerical value of the opacity is mainly 
dominated by the lack of precise knowledge of ^mo and Hq. For incoming neutrinos having 
the appropriate energy to interact resonantly, the opacity obtained ranges from ~ 0.05 to 
~ 0.3 for cosmological scenarios with 0.1 < ^2^/0 < 1 and cosmological constant parameter 
in the range < Qao < 0.7 [|l^. Although models in which the UHECR are produced by 



UHE neutrinos can require energy densities comparable to the luminosity of the Universe 
P), they would not necessarily have strong observable effects provided the opacity is low. If 
the opacities were larger, as one could expect for a mechanism generating the neutrinos at 
higher redshifts, there could be other observable consequences such as low energy photons 
above current experimental limits. Whatever the origin of the interacting neutrinos, if one 
requires a neutrino flux well exceeding that of other particles such as electrons, photons, 
protons and neutrons, models can be found which are consistent with the low energy photon 
flux bound as shown by the examples in Ref. 0. 

In order to explain the UHECR spectrum nearby the Earth the production rate within 
the absorption distance of the cosmic rays in the CMB (~ 50 Mpc) is fixed by data. If 
the local relic neutrino density is known this normalizes the neutrino flux. We leave a local 
density enhancement factor, 10^, to account for possible clustering effects and assume a 
halo radius D, otherwise the probability of interaction for a neutrino is very small and the 
neutrino flux needed to produce the cosmic rays must be enormous and in conflict with both 
energy considerations and experimental neutrino bounds. The survival probability for the 
incoming neutrino flux is given by a local opacity factor td^ integrating Eq.( [T|) to the halo 
limit taken to be D. This probability has a large resonance peak at E^'^'^ = M| /(2r7i;^) of 
width 6Ei, = EyV z/Mz where Mz and T z are the Z mass and width respectively. As long 
as D is below 50 Mpc the upper z limit is small z < 0.01 and the opacity at the resonant 
energy can be well approximated by the ("static") expression: 

TDu (Euo) ^ D a,p (2mpE,o) n^o ^ 1-3 10"^ 10« ^-j^. (2) 

The probability of interacting locally in the halo is given by 'Pi{Ey) = 1 — e""^^" ~ tou- As 
long as its size is not extremely large, D ^ 50 Mpc, the interaction probability, "P/, is small 
except for very large density enhancement factors. 

Taking the local neutrino flux entering the halo region to be (j){E^), the injection energy 
through resonant Z^ production is simply given by: 



£, = j dE,E,Vi{E,) (j){E, 



^ 1.3 [E-f vAE^r) 0(^r) ^, (3) 

where the last expression corresponds to the common approximation used for integrating the 
resonant cross section and the factor makes the expression numerically exact for a neutrino 
spectral index of 2. Following Waxman, Eq.(|^) can be equated to the produced energy flux 
of cosmic rays to obtain (pi^El^^): 

rEmax 

Sp= dEEME), (4) 

where (t>p{E) is the higher energy cosmic ray flux tail assumed to be due to this mechanism. 
For a neutrino spectral index 7, the flux is: 

' E,. 



m - ^' 



(5) 



1.3 [El^'Y 'PiiEl^') ^z/Mz [El^\ 
The important point is that 0(-E'^^*) is inversely proportional to the interaction probabil- 
ity at the resonance peak Vi{El^^) and to its width 6E1,. One should expect extrapolations 
of this flux with a fairly constant spectral index 7 both below and above the resonant energy 
because experience in astrophysical fluxes and theoretical considerations are very strong in 
supporting a neutrino flux spanning a few decades in energy. High energy fluxes from low 
interacting particles are severely constrained by existing experiments |[12|. For the range of 



energies considered here the strongest limit is given by the Fly's Eye group ||T9[. The non 



observation of horizontal air showers allows to put a limit on the integrated flux of any low 
interacting particle. Provided the neutrino flux can be extrapolated to the effective energy 
threshold for the Fly's Eye bound, Ep ~ 10^ GeV: 

^,{Ep) = / dEME) [1 - Vi{E,)] < $0. (6) 
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Fixing the neutrino mass, assuming 7 > 1, choosing a conservative (high) value of Emm = 
5 10^^ eV, and 2.5 for the proton spectral index in Eq.( ^, we can constrain the region of 
allowed values of 7 and Vi, using Eq.( |^) for (pi,. This is shown in fig. |l] for rrii, = 0.1, 1, 
and 10 eV. The figure shows that there is a critical spectral index 7 = 2.15 above which the 



model is ruled out for all masses in the 0.1 — 10 eV range. That is if 7 > 2.15 horizontal 
showers should have been observed even in the event that all neutrinos in the resonant 
energy range were converted to UHECR. If 7 ~ 1.2 however, a very low (depending on the 
neutrino mass) conversion probability could be allowed by data. 

Using Eq.( ^ to substitute the probability Vi{Ey) into the experimental limit expression 
in Eq. (P) we get a region of allowed parameter space ^, D for any given value of nii, and 7. 
This is shown in fig. § where the limits are given as the continuous lines for different values 
of the neutrino spectral index. 

Further restrictions apply in this parameter plot. The maximum density is constrained 



by the Fermi distribution to be [p,20 
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where v is the characteristic neutrino velocity in the halo. However, given the strong depen- 
dence on the neutrino mass and an unknown velocity we take it as a free parameter. 

In addition, if the total number of background neutrinos in the Universe is fixed, the 
density enhancement factors in the halo, their sizes and the maximum total number of halos 
are related. For constant density halos, assuming that no neutrinos are outside them, the 
number of neutrino clusters, in a Hubble radius, of a given size and a given enhancement 
factor is simply A^^ = 10~^(-R_f//-D)'^, where Rh is the Hubble radius. One can now easily 
see that the maximum number of halos can be read in a slant coordinate shown as dashed 
lines in Figs. ^J^. The three lines correspond to a single halo, 3 10^ and 10^° halos. The 
shaded region above the upper line implies less than 1 halo within a Hubble radius which is 
meaningless. The 3 10^ and lO^'' halo lines roughly correspond to one halo per supercluster 
(every 4 10^ Mpc'^) and to one per galaxy (every 500 Mpc^) respectively. Note that this is a 
maximum number of halos, for instance below the lower curve models can still have clusters 
around all galaxies, as long as the population of neutrinos in between the halos is non zero. 

Notice that protons are attenuated in the CMB in an energy loss distance of about 50 
Mpc. This means that for D > 50 Mpc the region of the halo outside a sphere of this radius 



centered around us can be ignored for the production of the local UHECR spectrum. Halo 
sizes exceeding 50 Mpc should be considered in these plots as having an effective size of 
50 Mpc. 

The approach is very conservative. Eq. (^) neglects fractions of the Z^ production energy 
which goes into particles that cannot be UHECR. The Z° decay will produce a particle flux 
following a typical fragmentation spectrum and the decay of the unstable particles will add 
low energy particles which cannot contribute to the UHECR. Neither can neutrinos from 
Z^ and pion decays nor that part of the high energy particles that are degraded by the 
showering developed in the intergalactic medium. The proton energy flux in Eq.(^ depends 
on the limits of the integration. It has been conservatively estimated by setting them close 
to the UHE part of the CR spectrum. As the observed cosmic ray flux spectral index at 
these energies is about 2.5, the lower integration limit gives the dominant contribution to 
the integral. The upper limit is not so relevant and it is in any case bounded by the neutrino 
resonance energy which in turn depends on the neutrino mass. For harder spectral indices 
closer to 2 one gets a similar result using an upper limit of order E^ax = 100 Emm as in 
Ref. @. 

We have also done an analytical calculation of the proton and photon fluxes originated 
from the neutrino-antineutrino annihilation again using Vi{E^) obtained from Eq. @). The 
flux of protons is given by: 

^^(E,^=J%^(e/-^^!^X, (8) 

where da/dEp is the cross section for the uu to produce a proton of energy Ep and X is 
the column depth of neutrinos in the halo. The cross section da/dEp can be written as the 
convolution: 

^ ^^ - / dy^— ^-^f{z/y), (9) 



dz Jz dy 

where z is the fraction of the energy taken by the proton, f{x) is the fragmentation function 
of a quark into an hadron, and da^uv -^ (l)/dy is the inclusive cross section for quark produc- 
tion from vv interactions. At the Z° resonance the non-resonant channels can be neglected at 
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this level of precision. We use Hill's fragmentation function ||2T|, f{x) = N 15/16x~^'^(l— x)^ 
with A^ = 0.03 for baryons. For calculating the secondary UHE photon spectrum similar 
expressions apply with A^ = 0.32 and an additional integral over vr*^ decay. These integrals 
are evaluated numerically. We normalize the photon plus proton spectrum to the observed 
cosmic ray flux at i^^ > 5 10^^ eV. We neglect the interaction of the high energy particles 
produced with the IR and CMBR which would increase the neutrino flux needed. 

If we again apply the Fly's Eye limit, parameter space becomes more restricted leaving 
less room for the conjecture as can be seen in fig. ^. No model with 7 > 2 is allowed 
in agreement with Ref. [§] but there is still room for harder spectral indices. The natural 
assumption of one halo per galaxy forcing small halo sizes and bounding the possible density 
enhancements is ruled out for any injection spectrum. This is unfortunate since a clear 
experimental signature of relatively small halo sizes would be a cosmic ray anisotropy due 
to the asymmetric position of the solar system within the halo, given by the ratio of our 
position and the galactic halo radius D. Assuming 10% sensitivity to anisotropy a future 
experiment such as the Auger Observatories |^ could test models with halos of order 



100 kpc. The picture is however not complete since neither Pauli blocking nor total mass 
constraints have been included. These must be related to mass density bounds that exist 
on different scales. 

The two case studies in Ref. use a halo size of 5 Mpc and density enhancements 
of 300 and 1000 (^ ~ 2.5 and ^ ~ 3). As stated in Ref. they are not in conflict with 
horizontal air shower data; this is because of the very hard spectral indices used, 7 ~ 1.2. 
It is interesting to notice, however, that the maximum number of halos corresponding to 
these models is ~ 7 10^ and ~ 2 10^ respectively, implying that most galaxies can not have 
an associated neutrino halo, yet the size of these halos is of order the average inter galactic 
distance. 

The total mass in neutrinos of such enhanced density halos is, on the other hand: 

M^ = -7rm^n^D^ = 1.2 10^° D^IO^ (~^) ^® (^°) 



Assuming a 1 eV neutrino mass this gives at least Mi, = 4.5 10^'^Mq for the most favourable 
case with density enhancement 300. Although this mass is suggestive of supercluster scales, 
its size fits closer to the Local Group. Unfortunately the case studies imply a total mass 
in neutrinos that clearly exceeds the dynamic mass measurements associated to the Local 
Group. For Q = 1 the mass within a sphere of radius R given by HqR = 390 km s~^ is 
determined to be Mq = 5.7 10^"^ Mq |2^. The mass increases by about 10% for Q = O.L 
Although the horizontal shower data and the halo mass constrains can be both met by a 
model neutrino mass (ttIj, < 0.1 eV) and spectral index 7 < 1.2, Pauli blocking arguments 
would imply orbital velocities exceeding 1000 km s~^ implying that a constant density model 
is inconsistent. 

The only other possibility that is left for these models is an even larger halo size. A 
halo size of ~ 50 Mpc, a neutrino spectral index of 7 ~ 1.2 and a neutrino mass of order 
0.1 eV seems viable. This case corresponds to a halo mass of order 10^^ Mq and orbiting 
velocities in the 1000 km~^ range. Increasing the halo size helps mostly by removing the 
strong constrain on the total mass. 

In summary Horizontal Shower limits, mass constraints and Pauli blocking mechanisms 
leave very little room for UHE neutrinos to explain the origin of the UHECR. Only very 
large halos, relatively low neutrino masses and unusually flat spectral indices are marginally 
allowed. Up to now we have assumed that there is absolute lepton number symmetry and 
in such cases the density parameter for the neutrinos is fixed by the neutrino mass: 

"" h^93eV ^ ^ 

It is however remarkable that if there is lepton number asymmetry, as recently suggested 
in Ref. [^, density enhancement comes rather naturally and is distributed uniformly over 
the whole Universe. Assuming that neutrinos are degenerate, Q^, ~ 0.01 and rrii, ~ 0.07 eV, 
they obtain an enhancement factor of ~ 30. Our analysis has been repeated considering the 
interactions within a sphere of ~ 50 Mpc, revealing that the model is completely consistent 
with the Fly's Eye limit on neutrino fluxes, provided the spectral index of the UHE neutrinos 
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satisfies 7 ^ 1.8. 

In any case there is no experimental signature provided by anisotropy. In the end however 
it is fortunate that these models can be further tested by experiment. A promising signature 
lies in the identification of photons as a significant component of the UHECR. This issue can 



be addressed by future experiments such as the Auger Observatories [^. Most importantly 
the fact that horizontal showers provide such a strict bound on these models also implies that 
future neutrino experiments, having much larger acceptance for neutrinos than Fly's Eye, 
should be able to detect the postulated UHE neutrino fluxes. Here the Auger Observatories 
may also play a role together with other high energy neutrino detectors in construction or 
planning stages. 
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FIG. 1. Excluded region in the Vj, 7 parameter space for three different values of the neutrino 
mass as marked. 




FIG. 2. Full lines mark the excluded regions by HS data (below the curves) in the D (Mpc), ^ 
parameter space for different values 7 assuming mi^= 1 eV. The neutrino flux has been normalized 
to the total energy in the UHECR (see text). The dashed lines represent the maximum number of 
halos provided the total number of relic neutrinos is fixed. The star marks the position of model 
in Ref.[9]. 




FIG. 3. Same as fig. 2, with neutrinos normalized to the UHECR flux (see text). 
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